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Abstract 

The relationship of Choptuik scaling to the scale invariance of Einstein's 

equation is explored. Ordinary dynamical systems often have limit cycles: 

periodic orbits that are the asymptotic limit of generic solutions. We show 

how to separate Einstein's equation into the dynamics of the overall scale and 

the dynamics of the "scale invariant" part of the metric. Periodicity of the 

scale invariant part implies periodic self-similarity of the spacetime. We also 

analyze a toy model that exhibits many of the features of Choptuik scaling. 
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I. INTRODUCTION 



Recently Choptuik has found scaling phenomena in gravitational collapse. He numer- 
ically evolves a one parameter family of initial data for a spherically symmetric scalar field 
coupled to gravity. Some of the data collapse to form black holes while others do not. There 
is a critical value of the parameter separating those data that form black holes from those 
that do not. The critical solution (the one corresponding to the critical parameter) has the 
property of periodic self similarity: after a certain amount of logarithmic time the profile 
of the scalar field repeats itself with its spatial scale shrunk. For parameters slightly above 
the critical parameter the mass of the black hole formed scales like (p — p*) 1 where p is 
the parameter, p* is its critical value and 7 is a universal scaling exponent that does not 
depend on which family is being evolved. Numerical simulations of the critical gravitational 
collapse of other types of spherically symmetric matter were subsequently performed. These 
include complex scalar fields, 0] perfect fluids, || axions and dilatons f|] and Yang-Mills 
fields. || In addition scaling has been found in the collapse of axisymmetric gravity waves. 
Thus scaling seems to be a generic feature of critical gravitational collapse. In some of 
these systems the critical solution has periodic self-similarity while in other systems it has 
exact self-similarity. 

These phenomena were discovered numerically, so one would like to have an analytic 
explanation for why systems that just barely undergo gravitational collapse behave in this 
way. Gundlach [ff| and Koike, Hara and Adachi [§] have explained the scaling of black hole 
mass analytically subject to the following assumptions: (i) the critical solution is periodically 
self-similar and (ii) the critical solution has exactly one unstable mode. This still leaves 
unexplained the mystery of why the critical solution is periodically self-similar. 

While periodic self-similarity is an unusual property for a dynamical system, periodicity 
is not. In fact, many dynamical systems have limit cycles, i.e. periodic trajectories that are 
approached asymptotically for a large class of initial conditions. What is it about Einstein's 
equation that gives rise to periodically self-similar solutions rather than periodic ones? One 
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feature of Einstein's equation is scale invariance: let g a b be a solution of the vacuum Einstein 
equation and let k be a positive constant. Then kg a ^ is a solution of the vacuum Einstein 
equation. The same property holds for the Einstein-scalar equation (the system studied by 
Choptuik). Let (g a i>, (ft) be a solution of the Einstein-scalar equation and let k be a positive 
constant. Then (kg a b, 0) is a solution of the Einstein-scalar equation. This feature of scale 
invariance suggests that in some sense the metric can be decomposed into an "overall scale" 
and a "scale invariant part" and that these two pieces have very different dynamics. In 
particular a periodically self-similar metric could be realized as a periodic scale invariant 
part of the metric. Also an exactly self-similar metric could be realized as a static scale 
invariant part of the metric. Thus if the dynamical system of the scale invariant part of the 
metric has a limit point i.e. a point in phase space that is approached asymptotically for a 
large class of initial conditions (a usual feature for dynamical systems) then the metric has 
a critical solution that is exactly self-similar. 

This program of finding scale invariant dynamics for the metric and matter fields has 
been carried out in the spherically symmetric case for various kinds of matter. In spherical 
symmetry the metric takes the form (T| 

ds 2 = -a 2 (r,t)dt 2 + a 2 {r,t)dr 2 + r 2 dVL 2 . (1) 

With t or r as the overall scale, the quantities a and a are scale invariant. The dynamics then 
becomes a set of differential equations for a, a and the scale invariant matter variables. These 
equations have been found for radiation fl, a massless scalar field 0, scalar electrodynamics 



and the Einstein- Yang-Mills system. 10 



Homogeneous cosmologies have also been treated using this sort of approach. Wainwright 



and Hsu [11] have found a set of scale invariant variables that describe the dynamics of 



homogeneous, anisotropic spacetimes. These variables have been useful both in the work of 



reference |TT] and in recent work of Rendall. [|T2] The connection between scale invariance 
and discrete self similarity has been noted in a study by Traschen of fluctuations about 
an extremal black hole. 



One would like to carry out this program of finding scale invariant metric variables in 
cases where there are no symmetries. The renormalization group framework of reference 
|| is sufficiently general to accomodate the case of scale invariance with no symmetries. 
However, the coordinate invariance of general relativity allows more than one possible set 
of scale invariant variables. One would like to find a set appropriate for critical graviational 



collapse. In the Euclidean approach to quantum gravity |L4J the metric is decomposed into 
a conformal factor and a determinant=l piece. However, in the Lorentzian case Einstein's 
equations are a dynamical system described by the ADM formalism. It therefore makes 
sense to have the overall scale determined by some dynamical condition. 

In section 2 we treat the dynamics of a toy model that has many of the features of 
Choptuik scaling. In section 3 we modify the ADM formalism to write Einstein's equation as 
a dynamical system of scale invariant quantities. Section 4 is a discussion of the implications 
of these results. 



II. TOY MODEL 



We now consider a toy model that has many of the features of Choptuik scaling. This 
model is constructed by adding an overall scale degree of freedom to a model known to have 
a limit cycle. The model is a dynamical system consisting of three functions of time: a, b 
and c. Define the quantity s by 
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where an overdot denotes derivative with respect to t. Choose the equations of motion of 
this system to be 
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Here e is a small positive constant. Note that this system has the property of scale invariance: 
if (a, b, c) is any solutions of the equations of motion and k is any positive constant then 
(ka, kb, kc) is also a solution of the equations of motion. This suggests that one might better 
understand this system by dividing its variables into an "overall scale" and a "scale invariant 
part." To make this idea precise define the quantities N, a and f3 by 

N = (abc) 1/3 , (6) 

a = \n(a/N) , (7) 

P = MWN) • (8) 

Note that a and (3 are scale invariant. Thus the quantity N is the overall scale of the system 
and (a,/?) is the scale invariant part of the system. Expressing equations (3-5) in terms of 
these new variables we find 

d ( N\ / 2 \ 

a + — \— J = — a + a + ea [1 — a ) , (9) 

* + s(£)=" + <i • (10) 
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— a — ft + — (^J^j = a — [3 + a — ea [\ — a 2 ^j , 
Adding equations (9-11) we find 

*U)= a - < 12 > 

Using this result in equations (9-10) we find 

a = — a + ea (l — a 2 ^j , (13) 

P = j3 . (14) 

Note that the system (a, /3) is a dynamical system independent of the overall scale N and 
with equations of motion given by equations (13-14). The general solution of equation (14) 



is 

P = \ (A) + vpa) z l + \ Wo ~ Vfso) e"* (15) 

where /3 = /3(0) and v@ = $(0). Note that as t — > oo we have /3 — >■ ±00 or depending 
on whether the quantity /3 + v@q is respectively positive, negative or zero. Now consider a 
generic one parameter family of initial data for the system (a, b, c) depending on parameter 
p. In general there will be a range of p for which the long time evolution of the system gives 
(3 — > 00, a range of p for which (3 —>■ — 00 and a critical value p* of p for which /3 — > 0. 
This is analogous to the case of Choptuik scaling where there is a range of p for which the 
late time evolution gives a black hole, a range of p for which the late time evolution gives 
flat space and a critical parameter p* separating these two ranges for which the late time 
evolution gives the Choptuik critical solution. 

Now consider the evolution equation (13) for a. This is the van der Pol equation. [15| 
It is well known that this equation has a stable limit cycle. A solution of equation (13) to 
zeroth order in e is 

a = 2 cos(t + O ) (16) 

where 0o is a constant. Furthermore any initial data sufficiently close to the trajectory given 
in equation (16) will approach that trajectory asymptotically as t — > 00. It then follows that 
for generic one parameter families of initial data for the (a,/3) system (sufficiently close 
to the limit cycle) there is a critical value of the parameter for which the trajectory is 
asymptotically periodic. 

Now consider the behavior, at late times, of the overall scale N in the critical solution. 
With a given by equation (16) we find using equation (12) 

N = n exp [ac< + 2 sin(t + <j> )] (17) 

where n and k are constants. That is, the overall scale is a periodic function multiplied by 
an exponential. The asymptotic behavior of the critical solution expressed in terms of the 
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variables (a, b, c) is then 



a = n exp 



Kt + 2V2 sin (t + 0o + vr/4) 



(18) 



6 = n exp [k£ + 2 sin (t + <f> )] , 



c = n exp 



/tt + 2v^ sin(t + 0o - tt/4) 



(19) 
(20) 



Thus the critical solution for (a, b, c) is periodically self similar. After a certain amount of 
time the solution repeats its behavior with only its overall scale changed. 



III. SEPARATION OF EINSTEIN'S EQUATION 

We now apply the ideas developed in the previous section to Einstein's equation. Start 
with the standard ADM formalism with, for simplicity, zero shift. The spacetime metric has 
the form 

ds 2 = -N 2 dt 2 + h ik dx i dx k . (21) 

Here h ik is the intrinsic metric of the t = const, slices. The lapse N can be chosen arbitrarily. 
The extrinsic curvature of the t = const, slices is 

K lk = ^d t h lk . (22) 

(Here d t denotes d/dt). We would like to keep the same time and space coordinates (t,x l ) 
when the spacetime metric is multiplied by an overall constant. Thus under the transfor- 
mation h ik — > kh ik for a constant k we want N — > \f~k~ N. An evolution equation consistent 
with this condition is 

d t N=\N 2 K . (23) 

Actually any constant could be chosen to replace the 1/3; but as we will see later the choice 
of 1/3 has other nice features. We now wish to replace h ik with a scale invariant quantity. 
Define 

h ik = N~ 2 h ik . (24) 
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This hik is the "scale invariant part" of the spatial metric. Using equations (22) and (23) 
we find 

d t h tk = 2K lk (25) 
where K ik , the "scale invariant part" of the extrinsic curvature, is defined by 

K lk = N~ l (K ik - l -Kh lk ^ . (26) 

Note that K ik is trace free. This is due to the presence of the factor 1/3 in equation (23). 
Thus the extrinsic curvature essentially has two parts: the trace and the scale invariant 
part. We now need to find an evolution equation for K ik and we would like to express this 
equation in terms of scale invariant quantities. One such quantity is 

Ui^N^diN . (27) 

The evolution equation for K ik is 

d t K ik = D i D k N + N (2K tp K k p — K K ik + R lk - ^R lk ) . (28) 

Here Di and ^Ri k are respectively covariant derivative and Ricci tensor of the metric h ik 
and R ik is the spacetime Ricci tensor. From equations (22), (23) and (28) straightforward 
but tedious algebra gives an evolution equation for K ik . Expressing that equation in terms 
of scale invariant quantities yields 

d t K ik = ~ 3 (NK) K ik + 2K ip K k p + R lk - ^R lk + 2D t u k - 2u^ k 

+ \ kk \ {Z) R + 2u p u p - 2D p u p - h pg R pg ] . (29) 

Here D-i and ^Ri k are respectively covariant derivative and Ricci tensor of the metric hi k 
and all indices are lowered and raised with hi k and its inverse. 

We now extract from Einstein's equation the equations of a dynamical system whose vari- 
ables are scale invariant quantities. For the vacuum case these variables are {h ikl K ik , uj k ). 
In the case where matter fields are present, and where the matter field equations are scale 
invariant, these metric variables must be supplemented with scale invariant matter fields. Is 



the set (hik, Kik, u>k), along with the matter fields, complete? That is, can the time deriva- 
tive of each variable be expressed in terms of the other variables? Clearly this holds for the 
time derivative of hik- m equation (29) the only questionable term is the one proportional 
to NK. However, we now show that the Hamiltonian constraint gives NK in terms of 
(hik, Kiki u>k) and the matter fields. The Hamiltonian constraint is 

®R + K 2 - K ik K ik = 2G F n"n v . (30) 

Expressing this constraint in terms of scale invariant quantities we have 

NK = ±\l 



K pq K pq + 2tu p u p + AD p u p - (3) R + 2h pq R pq - N 2 R] 1/2 . (31) 



It is also helpful to have an evolution equation for NK. From equations (22), (23) and (28) 
it follows that 

d t (NK) = - - (NK) 2 + 5D p uj p + Auj p uj p + h pq R m - {3) R . (32) 

In equations (31) and (32) all indices are lowered and raised with h ik and its inverse. We 
also need an evolution equation for oju- From equations (23) and (27) it follows that 

d t u k = ^b k (NK) . (33) 

Equations (25), (29) and (33) are the evolution equations for the dynamical system of scale 
invariant quantities (hik, K ik , u k ). Here the auxiliary quantity NK is given by the constraint 
equation (31) and its evolution is given by equation (32). 

Suppose that we have a solution for the dynamical system (hik, K ik , oj k ). What else 
is needed to determine the metric? Clearly the quantity N along with h ik is enough to 
determine the spacetime metric. However, the quantity oo k already contains most of the 
information about N. The only missing piece of information is the value of N at a single 
point of space as a function of time. Pick a spatial point x l and define N (t) = N(t, x l ). Then 
the quantity iVo along with a solution of the scale invariant dynamical system determines 
the spacetime metric. How does the quantity iVo evolve? Applying equation (23) at the 
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point x we find 



d t \nN = ±(NK)(x} ) ) 



(34) 



However, the quantity NK is determined by the dynamical system (h ik , Therefore, 



up to an overall constant scale (the value of N at some initial time t ) the spacetime metric 
is determined by the scale invariant dynamical system. 



follows that No is an exponential function multiplied by a periodic function. It then follows 
that the spacetime is periodically self-similar. If instead (ha-, Kik, w k ) are independent of 
time then iV is an exponential function of time and the spacetime is exactly self-similar. 

We now consider how to specify the scale invariant equations for one particular type of 
matter: the massless, minimally coupled scalar field. Here the Ricci tensor is 



Now suppose that the variables (h ik , K ik , uj k ) are periodic functions of time. Then it 



= 87rV M 0V^ 



(35) 



and the scalar field satisfies the wave equation 



V^W = 



(36) 



The field is scale invariant, so the full set of scale invariant quantities is (h ik , K ik , uj k , <fi). 



It follows that 



R ik = 87r Di<j)D k <j) 



(37) 



N 2 R = 8n [A0£>> - {d t <f)) 2 



(38) 



Expressed in terms of scale invariant quantities equation (36) becomes 




(39) 



IV. DISCUSSION 



The analysis of the previous section shows that the Einstein-scalar equation can be 
separated into equations for scale invariant quantities and an equation for the overall scale. 
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Periodic self similarity of the Choptuik critical solution is then equivalent to the existence 
of a limit cycle for the scale invariant system. Thus the somewhat odd property of periodic 
self similarity is reduced to the more familiar property of limit cycles of dynamical systems. 
(Correspondingly the property of exact self-similarity is reduced to the property of limit 
points of dynamical systems). 

What remains to be seen is whether the particular set of scale invariant variables used in 
section 3 is appropriate for the Einstein-scalar system. That is, we need to find out whether 
the variables (hik, Kn~, uik, <t>) are periodic functions of time in the Choptuik critical solution. 
This question is presently under study. If the answer is no, then one would need to search 
for a different set of scale invariant variables to describe critical gravitational collapse. 

One would also like to know what physical processes are responsible for the existence of 
limit cycles in the process of gravitational collapse. In the van der Pol equation, limit cycles 
come about due to energy dissipation (both positive and negative). The van der Pol system 
is essentially the harmonic oscillator with a small perturbation. For the harmonic oscillator 
energy E = (a 2 + « 2 )/2 one can show using equation (13) that on the average the energy 
increases if the amplitude of oscillations is below that of the limit cycle and decreases if the 
amplitude is above that of the limit cycle. Thus the process of energy dissipation drives 
arbitrary trajectories to the limit cycle. What is the analog of energy dissipation (both 
positive and negative) in the spherically symmetric Einstein-scalar system? When the scalar 
field is weakly gravitating it disperses at late times; so the energy in a fixed region tends 
to decrease. For a strongly self gravitating scalar field the self gravity tends to concentrate 
the energy of the field into ever smaller regions. It is these two competing effects that, 
depending on their relative strengths, combine to form field dispersion, black hole formation 
or the critical solution. What is needed is to find an "energy-like" quantity for the scale 
invariant system whose evolution corresponds to effects of concentration or dispersion of the 
field. 

What about the more general case of systems with axisymmetry or no symmetry at 
all? Here the equations become much more complicated and proving the existence of limit 
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cycles becomes much more difficult. Here too a promising approach is to find an energy 
like quantity that is conserved on the limit cycle. An examination of the behavior of such 
a quantity should yield new intuitions on how scaling arises. It is therefore likely that an 
examination of the scale invariant dynamical system will be a powerful tool in understanding 
Choptuik scaling. 
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